Dynamics of Bloch Oscillations in Disordered Lattice Potentials 
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We present a detailed analysis of the dynamics of Bloch oscillations of Bose-Einstein condensates in 
disordered lattice potentials. Due to the disorder and the interparticle interactions these oscillations 
undergo a dephasing, reflected in a damping of the center of mass oscillations, which should be 
observable under realistic experimental conditions. The interplay between interactions and disorder 
is far from trivial, ranging from an interaction-enhanced damping due to modulational instability for 
strong interactions, to an interaction-reduced damping due to a dynamical screening of the disorder 
potential. 



o 

(N 



5-H 



I 

O 

o 



> 
m 
rn 

O 

o 



X 



Bloch oscillations (BOs) constitute one of the most 
fundamental quantum phenomena for particles in peri- 
odic potentials. Under the influence of a constant force, 
particles in such potentials undergo an oscillatory motion 
instead of being linearly accelerated Although BOs 

are strongly hnked to the dynamics of electrons in solids, 
they have not been observed in bulk crystalhne materi- 
als so far, since lattice imperfections such as defects and 
phonons damp the coherent electronic motion before a 
single BO cycle is completed. The first observation of 
BOs was achieved in so-called semiconductor super lat- 
tices 0], which exhibit much larger periodicities. How- 
ever, disorder also leads to a fast decay of the BOs in 
these systems. 

On the contrary, ultra-cold gases in optical lattices 
provide perfect periodic potentials, with neither defects 
nor phonons. As a consequence, these systems open un- 
precented possibilities for the detailed analysis of quan- 
tum transport in lattices, and in particular for the ob- 
servation of BOs with very long lifetimes 0]. However, 
in spite of their perfect periodicity, optical lattices allow 
for the controlled introduction of different types of dis- 
order 0, @, 0, Si • Recent experiments have analyzed the 
effects of controlled disorder on the properties of ultra- 
cold gases [i, [13, ini, [il E Hi- These experiments 
have clearly shown that ultra-cold gases are indeed very 
promising systems for the analysis of the intriguing in- 
terplay between disorder and interparticle interactions. 

This paper is devoted to the analysis of this interplay in 
the BO dynamics of Bose-Einstein condensates (BECs) 
in tilted optical lattices. The dynamics of BOs is usually 
analyzed in terms of the Wannier-Stark (WS) energy lad- 
der in such systems. Any disorder introduces an unequal 
spacing in this ladder and thus leads to a damping of 
the BOs. In addition, the BOs are non-trivially modified 
by the interactions in the system. Strong interactions 
enforce damping, due to dynamical instability. How- 
ever, weak interactions can cause a dynamical screening 
of the disorder potential, prolonging the lifetime of the 
BOs Although quasi one-dimensional systems allow 
for a qualitative understanding of the physics involved. 



in typical experimental conditions radial excitations play 
a non-negligible role. In this paper we therefore first 
develop a qualitative understanding of the dynamics of 
BOs in the ID case and then extend our analysis to the 
experimentally relevant 3D case. Finally, we discuss a 
method for the observation of damped BOs by analyz- 
ing the momentum distribution in tinic-of-flight (TOE) 
measurements . 

Consider the quasi-lD case, in which a BEG at low 
temperature is so strongly confined by an harmonic trap 
of frequency uj± in the xy plane, that the chemical poten- 
tial fi is much smaller then the transverse level spacing 
/i << hLj±. Under these conditions, the ID dynamics of 
the condensate wave function <I> along the ^-axis, is given 
by the Gross-Pitaevskii-equation (GPE) [l^ 



2m dz'^ 



+ V{z) + g\<i>\^ 



(1) 



where m denotes the atomic mass, and g — 4:TTashuj± is 
the ID coupling constant, with the s-wave scattering 
length. V{z) denotes an external potential of the form 

V{z) = muj'^z'^/2 + sEr sin^ikz) -~Fz + Vdis{z). (2) 

i.e. a superposition of an axial harmonic trap with fre- 
quency w, a tilted potential with slope F, a disorder 
potential Vdis{z), and a lattice potential of periodicity 
d ~ T^/k, and depth s in units of the recoil energy 
Er = h^k^/2m. The condensate ground state in the su- 
perimposed harmonic and lattice potential serves as the 
initial state for our simulations of the dynamics. 

Figure [1] shows the averaged position {z{t)) = 
J(iz(f>*z$ for a Rubidium condensate, with u!±_ = 27r x 
200 Hz, d = 412.5 nm, s = 5 and Fd/Er = 0.05, for 
different particle numbers N, axial trap frequencies w, 
and disorder depths. As disorder potential, we consider 
Gaussian noise with correlation length L ~ 3.3 ^m. We 
define the disorder depth, Va as twice the standard de- 
viation from its mean value, according to [1, [lH. Fig- 
ure [1] demonstrates that even for situations where inter- 
actions are negligible (A^ = 1) disorder induces strong 
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BO damping. The role of the interactions is far from 
trivial, as shown in Fig. [Jl since they may enhance or 
significantly reduce the BO damping. We analyze this 
intriguing physics in detail below. 
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FIG. 1: Averaged position of a BEG undergoing BOs. Top 
frame: A*' = 350 particles, uj = 1 Hz, and disorder depths 
VA/Er = (dashed), 0.02 (solid), and 0.04 (dotted). Bottom 
frame: VA/Er = 0.02 for iV = 1, = 3.5 i/z (dotted), iV = 
350, uj = 7Hz (dashed) and N = 700, uj^WHz (solid). The 
trap frequencies were adjusted to match the initial wavepacket 
widths. 

The BO damping can be understood from a simplified 
analysis in the tight-binding regime. In this regime, we 
may project the condensate wavefunction into the Wan- 
nicr basis of localized states |n), where n labels the lattice 
sites. $„ = (n|<I>) is provided by the discrete non-linear 
Schrodinger equation (DNLSE)[l7| 



(3) 



where e„ denotes the external (possibly disordered) po- 
tential and U the on-site interaction energy. Hq is 
the Hamiltonian for the tilted lattice in the absence of 
both disorder and interactions, Ho^n = ^J^ni^^i+i + 
$ri-i) + Fdn^n, with hopping constant J, and tilting 
potential Fdn. In absence of disorder and interactions, 
the particles perform BOs with frequency lobo ^ Fd/h, 
and amplitude zbq = 2Jd/?iLL}BO- The cigenstatcs of Hq 



are the WS states [18[, \Wm) = J2n Jn-m\n), where J„ = 
Jn{zBo) is the Bessel function of first kind. The eigenen- 
ergies form the well-known WS-laddcr E'^ — mhujBO- 
We project the DNLSE in the WS basis Vm = (M^m|$) 



m' ,s,s' 



(4) 



and 4™ '™- 



where Bmm' = J2n^nJn- 

U J2n Jn-mJn-m' Jn-ni" Jn-Tu'" ■ If hu>BO IS mUCh larger 

than other energy scales we can neglect terms in ^ that 
introduce energy jumps larger than Tiloboi i-C we can 
employ rotating-wavc-approximation (RWA) arguments. 



In the RWA the disorder preserves pm = |^m|^, just pro- 
viding a shift E^ + Bmm- Note, that on a longer time 
scale of several BOs the RWA may fail and the disorder 
eventually leads to a transfer of population between WS 
states. The interactions can, even in the RWA, lead to 
a transfer of population between WS states. However, 
since pm ~ Pm whereas the phase ipm of tpm evolves as 
~ Pm , we may consider pm as being constant (at least 
at short time scales of few BOs) if the atomic wavepacket 
is sufficiently broad [2l|. In that case, the energy of the 
WS states becomes 

Em — rnflUJBO + -Bmm + 2 (r^Pm-2n + r^Pm-2n-l) , 



With ^ and 



(5) 

= ^n+i.n- Hence, the ener- 
gies of the WS states are not equidistantly spaced, the 
wavepacket undergoes dephasing, and as a consequence 
the BOs are damped. In addition, the BEG width ex- 
periences a breathing dynamics [2l|, as shown in Fig. [S] 
This figure shows (at least during the first BOs) a good 
agreement between our results from a direct simulation 
according to equa tion (jS ]) and those obtained using WS 
states ipmit) = \//Om(0)e~*-^"*/'', where Em is given by 




FIG. 2: Gondensate width in units of zbo obtained from 
direct simulation of ^ (dashed) and from the WS analysis 
discussed in the text (solid), for TujJbo = 0.08J, U/Tiujbo = 
3.7, and en/fiujBo = 0.1 cos 0.167m. 

Let us now discuss the intriguing role of the interac- 
tions on the BOs in more detail. In the lower panel of 
Fig. [U a stronger damping then in the single particle case 
can be observed for large nonlinearity (see the curve for 

= 700). This interaction-induced damping is related 
to the so-called dynamical instability [l^, [2o|. This in- 
stability occurs when the quasi-momentum reaches the 
outer parts of the Brillouin zone and small perturbations 
of the condensate wave function grow exponentially in 
time [lij . This mechanism becomes predominant with 
growing nonlinearity, strongly damping the BOs. 

On the contrary, the BO damping may be significantly 
reduced compared to the single particle case for weak 
nonlinearity (see the curve for A^ = 350). This effect 
is caused by an interaction-induced dynamical screening 
of the disorder [2^ and can be qualitatively better un- 
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derstood with an alternative semi-classical description 
of the BOs. In the regime of weak nonlinearity, we 
can assume that the dynamics occurs within the low- 
est Bloch band, and that the dynamical instability is ir- 
relevant on the time-scales considered. Let us denote 
the exact mean field potential obtained by solving the 
GPE (P) by Vmfiz,t) = and consider the ef- 

fective single particle problem, for a particle in the low- 
est Bloch band under the influence of Vmf{z, t), the tilt- 
ing force F and the disorder Vdis{z). The single par- 



ticle Hamiltonian is given by H^ff 



-2Jcos{kd) 



Fz + Vdisiz) + Vmfiz,t). In the quasi-momentum pic- 
ture, k k, and z = id/dk. Wc assume that the 
bare amplitude of the BOs, zbo is much smaller then 
the spatial spread of the initial wave function. This in 
turn implies a very narrow initial momentum distribu- 
tion, centered at fco, so that, in the absence of disor- 
der and nonlinearity (zero-order solution), z[t) ~ z(0) — 
Zbo cos[{Ft/h + ko)d\. The full Heisenberg equation for 
it reads dk/dt = F - {^Vd^s|^z){z) - {dV„^f / dz){z,t). 
We solve it pertubatively by inserting the zero order 
solution. Again, assuming a sharp initial momentum 
distribution we obtain z{t) ~ £(0) — z bo cos{^j^ + 

kod - dt'[{dVd^s/dz){z) - {dV„,f/dz)iz,t')]}. In 

order to calucalte the dephasing rate the latter ex- 
pression has to be averaged over the initial spread 
of z{0). A reasonable estimate of the rate is 7^ ~ 

dt'[{dVd.s/dz){z{o)) + {dv^f/dz){m,t'm- 

Note that, when acting alone, both disorder and nonlin- 
earity lead to the damping of the BOs. However, when 
acting together, they may compensate each other if the 
product of the time averaged forces due to disorder and 
nonlinearity averaged over z(0) is negative, qualitatively 
explaining the dynamical screening of the disorder ob- 
served in Fig. [TJ 
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FIG. 3: Expectation value of the axial velocity for disorder 
depths (dashed), 0.02 (solid) and 0.06 (dotted), for 
the 3D example discussed in the text. 



In the previous discussion we have constrained our 
analysis to the somewhat simplified quasi-lD regime, 
where the xy dynamics is absent. However, typical exper- 



iments are not performed in this regime, and transversal 
excitations of the condensates may signficantly alter the 
BO dynamics. Hence a quantitative description of the 
BOs demands a three-dimensional GPE simulation. We 
maintain a cylindrical trap, and hence we simplify our 
3D calculations by assuming cylindrical symmetry of the 
wavefunction around the z-axis. Since experiments typ- 
ically detect BOs by observing the velocity distribution, 
we analyze the expectation value of the axial velocity 
< Vz >= J d^r^*-fi^dz^. Figure [3] shows the case of 
iV = 5 X IC^ particles, for the same lattice tilting con- 
sidered before, but for trap frequencies ojz — in x 14Hz, 
LUj_ = 2tt X 35Hz, lattice depth s ~ 2, and a disorder 
potential of depth Va = 0.02 and correlation length 
L lO^m. These parameters are well attainable ex- 



perimentally 111, and hence our simulations show that 



the damping of BOs can be studied under realistic con- 
ditions. For the large number of atoms considered here, 
damped BOs are observed even in the absence of dis- 
order due to a fast damping by the dynamical instabil- 
ity. Under these conditions our 3D results clearly deviate 
from the expected quasi-lD results, due to the develop- 
ment of a complex radial dynamics shown in Fig. [4l due 
to instability of the radial excitations. Hence, although 
the dynamical instability is also encountered in ID, the 
full computation of the dynamics does require the careful 
consideration of the radial degree of freedom. 
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FIG. 4: Radial density profiles at the axial position zo of 
maximal density, aX t = (solid), 16ms (dotted) for the case 
discussed in the text. 

Finally, we discuss a method for the experimental ob- 
servation of damped BOs. A direct measurement of the 
damping in position space constitutes a very difficult 
task, since the oscillation amplitude is typically too small 
for an in-situ measurement. TOF measurements provide 
a much better method. Due to the very fast decrease 
of the mean-field energy during the TOF, the expanded 
density distribution provides an approximative image of 
the momentum distribution n(p) of the sample at the mo- 
ment of release. Figure [5] shows the momentum distribu- 
tion at various stages during the BOs. In the absence of 
disorder, the quasi-momentum q scans the Brillouin-zone 
due to the acceleration introduced by the tilting force 
(q — F/h) [1^. Consequently, the population of the dif- 
ferent momentum components changes, resulting in a co- 
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FIG. 5: Gray-scale plot of the momentum distribution n{p± — 
0,Pz), calculated every 1 ms during the BOs. The upper panel 
shows the non-disordered case, the lower panel shows the re- 
sult for a disorder depth of 0.06 -Er- The pictures correspond 
to the 3D configuration of Figs. [3] and U) 



herent oscillation of the mean momentum. The spectrum 
consists of several sharp peaks, which are separated from 
each other by the lattice momentum ps = ^hk. Hence, 
the density distribution after an expansion time r con- 
sists of several small clouds, well separated from each 
other by Az ~ 2pBT/m. This picture changes when dis- 
order is introduced to the system. The initial sharp mo- 
mentum components are progressively broadened, even- 
tually reaching an irregular occupation of momenta. This 
broadening originates from the irregular energy spacing 
of the WS states and leads to a broadened density dis- 
tribution after TOF. Note, that on long time scales, a 
significant broadening of the momentum components is 
also introduced by the interactions in the absence of dis- 
order, as shown in Fig. [Sj However, a very clear distinc- 
tion between the ordered and disorderd case is possible 
on time scales of up to a few BO periods. Moreover, 
our 3D simulations of the expansion dynamics show that 
it is possible to observe a clear separation between the 
center of mass position for the ordered and disordered 
lattices under realistic conditions, due to the differences 
in the expected value of the momentum in both cases 
(Fig. [3]). Also note that Fig. [5] shows the axial momen- 
tum distribution for zero transversal momentum, which 
is significantly reduced during the time evolution. This 
reflects the previously mentioned excitation of transverse 
modes. 

In summary, both disorder and interactions separately 



lead to BO damping. However, when acting simultane- 
ously the interactions may partially screen the disorder, 
leading to a reduction of the BO damping. We have 
shown that the interplay of disorder and interactions may 
be observed under realistic experimental conditions by 
monitoring the evolution of the momentum distribution 
of the system. Although we were mainly interested in dis- 
ordered lattices, similar results may be obtained for lat- 
tices in the presence of inhomogeneous forces, as e.g. spa- 
tially inhomogeneous Casimir-Polder forces close to sur- 
faces. It has recently been proposed that the frequency 
shift of BOs of lattice fermions can be an excellent way 
of measuring such tiny forces (23j . The monitoring of BO 
damping and the sample width may provide an excellent 
way of proving the inhomogeneity of these forces. 
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